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Abstract 



O \ The probability of a charged particle production by the electric field of a charged 

black hole depends essentially on the particle energy. This probability is found in the 
nonrelativistic and ultrarelativistic limits. The range of values for the mass and charge 
of a black hole is indicated where the discussed mechanism of radiation is dominating 
over the Hawking one. 
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1. The problem of particle production by the electric field of a black hole has been 
discussed repeatedly [1-6]. The probability of this process was estimated in these papers 
using in some way or another the result obtained previously [7-9] for the case of an electric 
field constant all over the space. This approximation might look quite natural with regard 
to sufficiently large black holes, for which the gravitational radius exceeds essentially the 
Compton wave length of the particle A = 1/m. (We use in the present paper the units with 
h = 1, c — 1; the Newton gravitational constant k is written down explicitly.) However, in 
fact, as will be demonstrated below, the constant-field approximation, generally speaking, is 
inadequate to the present problem, does not reflect a number of its essential peculiarities. 

It is convenient to start the discussion just from the problem of particle creation by a 
constant electric field. Here and below we restrict to the consideration of the production of 
electrons and positrons, first of all because the probability of emitting these lightest charged 
particles is the maximum one. Besides, the picture of the Dirac sea allows one in the case of 
fermions to manage without the second-quantization formalism, thus making the considera- 
tion most transparent. 

To calculate the main, exponential dependence of the effect, it is sufficient to restrict to a 
simple approach due to |7j (see also the textbook JjlJ). In the potential —eEz of a constant 
electric field E the usual Dirac gap (Fig. 1) tilts (see Fig. 2). As a result, a particle which 
had a negative energy in the absence of the field, can now tunnel through the gap (see the 
dashed line in Fig. 2) and go to infinity as a usual particle. The hole created in this way is 
nothing but antiparticle. An elementary calculation leads to the well-known result for the 
probability of particle creation: 



First of all, the action inside the barrier does not change under a shift of the dashed line 
in Fig. 2 up or down. Just due to it expression (|lj) is independent of the energy of created 
particles. Then, for the external field to be considered as a constant one, it should change 
weakly along the path inside the barrier. However, the length of this path is not directly 
related to the Compton wave length of the particle. In particular, for an arbitrary weak field 
the path inside the barrier becomes arbitrary long. 

Thus, one may expect that the constant-field approximation is not, generally speaking, 
applicable to the problem of a charged black hole radiation, and that the probability of 
particle production in this problem is strongly energy-dependent. The explicit form of this 
dependence will be found below. We restrict in the present work to the case of a non-rotating 
black hole. 

2. We start the solution of the problem with calculating the action inside the barrier. 
The metric of a charged black hole is well-known: 




(1) 



This simple derivation explains clearly some important properties of the phenomenon. 



ds 2 = fdt 2 



f- 1 dr 2 -r 2 {d9 2 + sin 2 0# 2 ), 
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M and Q being the mass and charge of the black hole, respectively. The equation for a 
particle 4-momentum in these coordinates is 



r l («-f)' -/!>"-£=.»■. w 

Here e and p are the energy and radial momentum of the particle. We assume that the 
particle charge e is of the same sign as the charge of the hole Q, ascribing the charge — e to 
the antiparticle. 

Clearly, the action inside the barrier is minimum for the vanishing orbital angular momen- 
tum I. It is rather evident therefore (and will be demonstrated in the next section explicitly) 
that after the summation over I just the s-state defines the exponential in the total proba- 
bility of the process. So, we restrict for the moment to the case of a purely radial motion. 
The equation for the Dirac gap for I = is 

e±(r) = —±mJf. (5) 
r v 



It is presented in Fig. 3. It is known 1] that at the horizon of a black hole, for r = r + = 
kM + ^Jk 2 M 2 — kQ 2 , the gap vanishes. Then, with the increase of r the lower boundary 
of the gap e_(r) decreases monotonically, tending asymptotically to —m. The upper branch 
e + (r) at first, in general, increases, and then decreases, tending asymptotically to m. 

It is clear from Fig. 3 that those particles of the Dirac see whose coordinate r exceeds 
the gravitational radius r + and whose energy e belongs to the interval e_ (r) > e > m, tunnel 
through the gap to infinity. In other words, a black hole looses its charge due to the discussed 
effect, by emitting particles with the same sign of the charge e, as the sign of Q. Clearly, the 
phenomenon takes place only under the condition 

e< 2 ^ ffl \ 
— > m. (o) 

r + 

For an extreme black hole, with Q 2 = kM 2 , the Dirac gap looks somewhat different (see 
Fig. 4): when Q 2 tends to kM 2 the location of the maximum of the curve e+(r) tends to r + , 
and the value of the maximum tends to eQ /r + . It is obvious however that the situation does 
not change qualitatively due to it. Thus, though an extreme black hole has zero Hawking 
temperature and, correspondingly, gives no thermal radiation, it still creates charged particles 
due to the discussed effect. 

In the general case Q 2 < kM 2 the doubled action inside the barrier entering the expo- 
nential for the radiation probability is 

rrn 

2S = 2 dr\p(r,e)\ 
J n 

= 2 jT r2 _ 2k d Z + kQ 2 + 2(eeg ~ km2M)T ~ (e2 ~ km2)Q2 - (?) 
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Here po = \/e 2 — m 2 is the momentum of the emitted particle at infinity, and the turning 
points 7-1,2 are as usual the roots of the quadratic polynomial under the radical; we are inter- 
ested in the energy interval m < e < eQ/r + . Of course, the integral can be found explicitly, 
though it demands somewhat tedious calculations. However, the result is sufficiently simple: 



2 

m 



2S = 2tt— — [ e Q-(e-p )kM\. (8) 

(e + Po) Po 

Certainly, this expression, as distinct from the exponent in formula ([I]), depends on the energy 
quite essentially. 

Let us note that the action inside the barrier does not vanish even for the limiting value 
of the energy e m = eQ/r + . For a nonextreme black hole it is clear already from Fig. 3. For 
an extreme black hole this fact is not as obvious. However, due to the singularity of \p(r, e)\, 
the action inside the barrier is finite for e = e m = eQ/r + for an extreme black hole as well. 
In this case the exponential factor in the probability is 

exp |^-7r kmM^j . (9) 

Due to the extreme smallness of the ratio 

\fkm 



10~ 21 , (10) 

the exponent here is large only for a very heavy black hole, with a mass M exceeding that of 
the Sun by more than 5 orders of magnitude. And since the total probability, integrated over 
energy, is dominated by the energy region e ~ e m , the semiclassical approach is applicable in 
the case of extreme black holes only for these very heavy objects. Let us note also that for 
the particles emitted by an extreme black hole, the typical values of the ratio e/m are very 
large: 

e Cm e Q 6 1 n 21 

— ~ = = — j= — ~ 1U . 

m m kmM ykm 
In other words, an extreme black hole in any case radiates highly ultrarelativistic particles 
mainly. 

Let us come back to nonextreme holes. In the nonrelativistic limit, when eQ/r + — > m 
and, correspondingly, the particle velocity v — > 0, the exponential is of course very small: 

exp I I. (11) 

Therefore, we will consider mainly the opposite, ultrarelativistic limit where the exponential 
is 

-7i — eQ\. (12) 

Of course, here also the energies £ ~ € m ~ cQ/kAd are essential, so that the ultrarelativistic 
limit corresponds to the condition 

eQ > kmM. (13) 



3 



But then the semiclassical result ( |I2"D is applicable (i.e., the action inside the barrier is large) 
only under the condition 

kmM > 1. (14) 

Let us note that this last condition means that the gravitational radius of the black hole 
(r_|_ ~ kM) is much larger than the Compton wave length of the electron 1/m. In other 
words, the result ( |T2"D refers to macroscopic black holes. Combining ( |T3"D with (|T4|), we arrive 



at one more condition for the applicability of formula (|T2|): 



eQ > 1. (15) 

We will come back to this relationship below. 

Let us note that in Q the action inside the barrier was being calculated under the same 
assumptions as formula (|T2|). However, the answer presented in [|J, 25 = 7rm 2 r 2 h /eQ, is 
independent of energy at all (and corresponds to formula (|I|) which refers to the case of a 
constant electric field). I do not understand how such an answer could be obtained for the 
discussed integral in the general case e^e m . 

3. The obtained exponential is the probability that a particle approaching the turning 
point 7™i (see Figs. 3, 4) from the left, will tunnel through the potential barrier. One should 
recall that in the general case the position of the turning point depends not only on the 
particle energy e, but on its orbital angular momentum I as well. The total number of 
particles with given e and /, approaching a spherical surface of the radius r% in unit time, is 
equal to the product of the area of this surface 

S = Att rl{e, I) (16) 

times the current density of the particles 



/Pbo dt 

(see, e.g., W2\, §90). The particle velocity is as usual 

'-£-1 < l8 > 

(the subscript r of the radial momentum p is again omitted). To obtain an explicit expression 
for the particle density p, we will use the semiclassical approximation (the conditions of its 
applicability for the region r + < r < n will be discussed later). Let us note that the volume 
element of the phase space 

(19) 



dp x dp y dp z dxdydz 



(2tt) 3 

is a scalar. (The factor 2 here is due as usual to two possible orientations of the electron 
spin.) On the other hand, the number of particles in the elementary cell dxdydz equals (see 
0,§9O) 

Py/^dxdydz, (20) 
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where 7 is the determinant of the space metric tensor. Since all the states of the Dirac sea 
are occupied, we obtain by comparing formulae ([19]) and ( f20|) that the following expression 

p 2 ^ dp x dp y dp z _ 2 ^ dp x dp y dp z 



9oo V#oo7 (2vr) 3 yf=g ^ (2tt) 



should be plugged in formula ( |T7| ) for the current density (the summation here and below is 
performed with fixed e and I, see (|17D). In our case the determinant g of the four-dimensional 
metric tensor does not differ from the flat one, so that the radial current density of the 
particles of the Dirac sea is 

The summation in the right-hand-side reduces in fact to the multiplication by the number 
21 + 1 of possible projections of the orbital angular momentum 1 onto the z axis and to the 
integration over the azimuth angle of the vector 1, which gives 2n. With the account for the 
identity 

-^—dp r = de, 
dp r 

we obtain in the result 

fie I) = 2 M21 + 1} (22) 

Finally, the pre-exponential factor in the probability, differential in energy and orbital angular 
momentum, is 

2(2/ + 1) 



7T 

Correspondingly, the number of particles emitted per unit time is 



(23) 



dN 2 
dt 71 



yde^(2Z + l)exp[-2S(e,Z)]. (24) 



In the most interesting, ultrarelativistic case dN/dt can be calculated explicitly Let us 
consider the expression for the momentum in the region inside the barrier for / 7^ 



\p(e,l,r)\ = f- 1 



\ 



The main contribution to the integral over energies in formula ( p4]) is given by the region 
e — > e m . In this region the functions /(r) and e — eQ/r, entering expression (p5|) , are small 
and change rapidly. As to the quantity 

/i 2 (r,/)=m 2 + ^, (26) 
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one can substitute in it for r its average value, which lies between the turning points r% and 
r 2 . Obviously, in the discussed limit e — > e m the near turning point coincides with the horizon 
radius, r\ — r + . And the expression for the distant turning point is in this limit 



r 2 = r + 



2/i 2 y/k 2 M 2 - kQ 2 



(27) 



Assuming that for estimates one can put in formula ( PB| ) r ~ r + , one can easily show that 
the correction to 1 in the square bracket is bounded by the ratio l 2 /(eQ) 2 . Assuming that 
this ratio is small (we will see below that this assumption is self-consistent), we arrive at 
the conclusion that r 2 « r+, and hence /i 2 can be considered independent of r: /x 2 (r, /) = 
m 2 + I 2 jr\. As a result, we obtain 




2%,0~7reQ + . (28) 



Now we find easily 



oW ( eQ \ 
— — = m 



3 / 2 2 

7 urn r 



, exp -± . (29) 

at v 7rmr + y \ eQ J 

Let us note that the range of orbital angular momenta, contributing to the total probability 
(p5D, is effectively bounded by the condition I 2 < eQ. Since eQ ^> 1, this condition allows 
one to change from the summation over I in formula ( p4]) to the integration. On the other 
hand, this condition justifies the used approximation fi 2 (r, I) = m 2 + I 2 jr\. 

However, up to now we have not considered one more condition necessary for the deriva- 
tion of formula (|29|). We mean the applicability of the semiclassical approximation to the 
left of the barrier, for r + < r < t\. This condition has the usual form 

dr p(r) ^ 

In other words, the minimum size of the initial wave packet should not exceed the distance 
from the horizon to the turning point. Using the estimate 

r + (eQ - er + ) 
[r — r + )(r — r_J 

for the momentum in the most essential region, one can check that for an extreme black hole 
the condition ( p0|) is valid due to the bound eQ ^> 1. In a non-extreme case, for r + — r_ ~ r + , 
the situation is different: the condition ( ^0|) reduces to 

eQ - 1 eQ 
e < ~ — . (31) 

r + r + 

Thus, for a non-extreme black hole in the most essential region e — > e m the condition of the 
semiclassical approximation is not valid. Nevertheless, the semiclassical result ( |2H| ) remains 
true qualitatively, up to a numerical factor in the pre-exponential. 



6 



In concluding this section few words on the radiation of light charged black holes, for 
which kmM < 1, i.e., for which the gravitational radius is less than the Compton wave 
length of the electron. In this case the first part, 

eg - 1 

e < , 

r+ 

of inequality ([H]), which guarantees the localization of the initial wave packet in the region 
of a strong field, means in particular that 

eQ = Za > 1 (32) 



(we have introduced here Z = Q/e). It is well-known (see, e.g., JL3], jlUfl ) that the vacuum 
for a point-like charge with Za > 1 is unstable, so that such an object looses its charge 
by emitting charged particles. It is quite natural that for a black hole whose gravitational 
radius is smaller than the Compton wave length of the electron, the condition of emitting 
a charge is the same as in the pure quantum electrodynamics. (Let us note that the unity 
in all these conditions should not be taken too literally: even in quantum electrodynamics, 
where the instability condition for the vacuum of particles of spin 1/2 is for a point-like 
nucleus just Za > 1, for a finite-size nucleus it changes fll3| , |14| to Za > 1.24. On the other 
hand, for the vacuum of scalar particles in the field of a point-like nucleus the instability 



condition is jrj, |TjJ: Za > 1/2.) As has been mentioned already, for a light black hole, 
with kmM < 1, the discussed condition eQ > 1 leads to a small action inside the barrier 
and to the inapplicability of the semiclassical approximation used in the present article. The 
problem of the radiation of a charged black hole with kmM < 1 was investigated numerically 

m 



in 



4. The exponential 

exp 



2 2 

■Km ri 



eQ J 

in our formula ( p9|) coincides with the expression arising from formula ([!]), which refers to a 
constant electric field E, if one plugs in for this field its value Q/r 2 + at the black hole horizon. 
As has been mentioned already, an approach based on formulae for a constant electric field 
was used previously in Refs. [1-6]. Thus, our result for the main, exponential dependence 
of the probability integrated over energies, coincides with the corresponding result of these 
papers. Moreover, our final formula (|4]) agrees with the corresponding result of Ref. |§ up to 
an overall factor 1/2. (This difference is of no interest by itself: as has been noted above, for 
a non-extreme black hole the semiclassical approximation cannot guarantee at all an exact 
value of the overall numerical factor.) 

Nevertheless, we believe that the analysis of the phenomenon performed in the present 
work, which demonstrates its essential distinctions from the particle production by a constant 
external field, is useful. First of all, it follows from this analysis that the probability of the 
particle production by a charged black hole has absolutely nontrivial energy spectrum. Then, 
in no way are real particles produced by a charged black hole all over the whole space: for 
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a given energy e they are radiated by a spherical surface of the radius r2(e), this surface 
being close to the horizon for the maximum energy. (It follows from this, for instance, that 
the derivation of the mentioned result of Ref. M for dN/dt has no physical grounds: this 
derivation reduces to plugging E = Q/r 2 into the well-known Schwinger formula ||, obtained 
for a constant field, with subsequent integrating all over the space outside the horizon.) 

Let us compare now the radiation intensity / due to the effect discussed, with the intensity 
Ih of the Hawking thermal radiation. Introducing additional weight e in the integrand of 
formula (E4T), we obtain 



/= w(^Vexp(-^). (33) 
\ixmr + J \ eQ J 

As to the Hawking intensity, the simplest way to estimate it, is to use dimensional arguments, 
just to divide the Hawking temperature 

1 



by a typical classical time of the problem r + (in our units c = 1). Thus, 

1 



2 ■ 



(34) 



More accurate answer for Ijj differs from this estimate by a small numerical factor ~ 2 • 10 -2 , 
but for qualitative estimates one can neglect this distinction. The intensities (|33p and fl3~4] ) 
get equal for 

eQ ~ - (mr+)2 ~ 1 (fcmM)2 . (35) 
6 ln(mr + ) 6 \n(kmM) 

(One cannot agree with the condition eQ ~ l/(47r) for the equality of these intensities, 
derived in Ref. from the comparison of e m = eQ/r + with T H = l/(47rr + ) .) 

Let us consider in conclusion the change of the horizon surface of a black hole, and hence 
of its entropy, due to the discussed non-thermal radiation. To this end, it is convenient to 



introduce, following Ref. fHf] , the so-called irreducible mass M of a black hole: 



2M = M + \jM 2 - Q 2 ; (36) 
here and below we put k — 1. This relationship can be conveniently rewritten also as 

M = M " + & a < 37 > 

Obviously, r + = 2Mo, so that the horizon surface and the black hole entropy are proportional 
to M 2 . 

When a charged particle is emitted, the charge of a black hole changes by AQ = — e, and 
its mass by AM = —eQ/r + + £, where £ is the deviation of the particle energy from the 
maximum one. Using the relationship fl37D , one can easily see that as a result of the radiation, 
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the irreducible mass Mo, and hence the horizon surface and entropy of a non-extreme black 
hole do not change if the particle energy is the maximum one eQ/r + . In other words, such 
a process, which is the most probable one, is adiabatic. For £ > 0, the irreducible mass, 
horizon surface, and entropy increase. 

As usual, an extreme black hole, with M = Q = 2Mo, is a special case. Here for the 
maximum energy of an emitted particle e m = e, we have AM = AQ = — e, so that the black 
hole remains extreme after the radiation. In this case AMq = — e/2, the irreducible mass 
and the horizon surface decrease. In a more general case, AM = — e + £, the irreducible mass 
changes as follows: 



Clearly, in the case of an extreme black hole of a large mass, already for a small deviation £ of 
the emitted energy from the maximum one, the square root is dominating in this expression, 
so that the horizon surface increases. 
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